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Abstract. We show that high quality laser guide star (LGS) adaptive optics (AO) observations of 
nearby early-type galaxies are possible when the tip-tilt correction is done by guiding on nuclei 
while the focus compensation due to the changing distance to the sodium layer is made 'open loop'. 
We achieve corrections such that 40% of flux comes from R<0.2 arcsec. To measure a black hole 
mass (M,) one needs integral field observations of both high spatial resolution and large field of 
view. With these data it is possible to determine the lower limit to M, even if the spatial resolution 
of the observations are up to a few times larger than the sphere of influence of the black hole. 
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INTRODUCTION 

The correlation between the masses of supermassive black holes (SMBH) with the 
global characteristics of their host galaxies [e.g. 1, 2] is crucial for our interpretation 
of galaxy formation and evolution [e.g. 3, 4]. M, is usually determined using dynamical 
tracers such as a disk of gas clouds in Keplerian rotation around the SMBH or the line- 
of-sight velocity distribution (LOSVD) of stars. Both methods require observations of 
high spatial resolution since the dynamical models have to be constrained by resolved 
kinematics which probe the sphere of gravitational influence of the SMBH (Rsph), 
defined as the distance from SMBH at which the potential of the galaxy and SMBH 
are approximately equal. This spatial scale, usually defined as Rsph = GM,/o^, where 

is the velocity dispersion of the galaxy, is typically significantly less than an arcsec in 
apparent size, even for nearby galaxies. 

The influence of the SMBH, however, will be felt at larger radii than Rshp, albeit to a 
lessening degree. Here we show that it is possible to constrain the M, when Rshp is up to 
3 times smaller than the spatial resolution of the observations if the observation consists 
of (i) low (arcsec) spatial resolution integral field data covering the galaxy out to about 

1 effective radius required to determine the overall dynamical mass-to-light ratio of the 
system, ( ii) high (sub-arcsec) resolution integral field data probing the stellar kinematics 
in the vicinity of the SMBH and have (Hi) sufficient signal-to-noise ratio (S/N) to extract 
higher order moments of the LOSVD. 



In this report we summarise the results of [5] which present the usage of a new 
observational method with LGS AO system especially suitable for determination of M, 
in the nuclei of nearby low mass early-type galaxies. For determination of M, in massive 
early-type galaxies using similar method, but non-AO observations see the contribution 
of [6] in this proceedings. 

MEASURING SMALL M. FROM GROUND 

The main prohibitive issue with AO observations is that, even with LGS capabilities, a 
natural guide star is required for an optimal correction of atmospheric aberrations and 
there are only a handful of galaxies with a near-resolvable Rsph that have a suitable guide 
star. This limitation has prompted innovative use of AO techniques, such as neglecting 
altogether the low-order corrections provided by the natural tip-tilt [7]. In this study, we 
employ a novel method developed at Gemini Observatory for Altair LGS AO system. 



Open loop focus model of LGS 

To avoid the issue of suitable guide stars, it is, in principle, also possible to use the 
nucleus of the galaxy as a natural guide source for the system, provided it is sufficiently 
bright and compact (drop of > 1 magnitude within the central ~ 1 arcsec for Altair). The 
central light profiles of early-type galaxies show a general change from steeply rising 
'cusp' profiles at lower masses, towards a flatter central profile that can define a central 
'core' region [e.g. 8]. In order to test what correction of the PSF could be expected for 
typical early-type galaxies, we chose NGC524 and NGC2549 which both satisfy the 
basic Altair requirement for nuclear guiding. Specifically, NGC524 has a core-like light 
profile, while NGC2549 a cusp-like profile at the HST resolution. In this proceeding we 
will consider the observations of NGC2549 only, the smaller of the two galaxies. 

After initial observations were attempted, it became clear that the chosen nuclei, 
while suitable for tip-tilt correction, are too faint to be used for constraining the focus. 
Gemini staff implemented a procedure by which the focus correction during the science 
integration is controlled by a geometric function that takes into account the change in the 
distance to the sodium layer as the telescope position changes. In this, so-called 'open- 
loop' focus model, the only time-dependent parameter is the altitude of the sodium layer, 
which is determined immediately before observing the galaxy by 'tuning' the LGS AO 
system using a nearby bright star. When all the control loops of the system (tip-tilt, 
focus, and LGS) have converged with this reference source, the loops are opened and 
the science target is acquired. The tip-tilt and LGS control loops are then closed, using 
the galaxy nucleus and laser beacon respectively as reference sources. The focus loop 
is left open, being passively controlled by the open-loop model. After approximately 
one hour of science observations, the bright reference star is re-observed, so that the 
degradation of the PSF can be estimated (our tests show minimal degradation) and the 
LGS AO system can be re-optimised for further observations. 

We estimated the PSF of the science observations by convolving (with a double 
Gaussian) and rebining (to the pixel size of our observations) an HSTAVFPC2 image. 
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FIGURE 1. Schwarzschild dynamical models and the determination of the best fitting parameters for 
NGC2549, Each symbol is a dynamical model. The agreement between the data and the models are 
described by overploted contours Ax^ contours showing 1, 2 and 3 a levels for two parameters. Further 
contours are spaced by a factor of 2, Large symbol on each panel marks the best fitting models. From 
left to right panels show grids of models constrained by both SAURON and NIFS kinematics, models 
constrained by SAURON only kinematics and models constrained only using NIFS data. 
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FIGURE 2. Comparison between data and model velocity dispersion maps for NGC2549. From left 
to right: a (symmetrised) and model prediction maps for different M, (at the best fitting M/L), The best 
fitting model is shown in the map adjacent to the observed a map, followed by a model with a too small 
M, and a model with too large M,, The models were constrained using both SAURON and NIFS data. 
Note the colour change in the central few pixies which are a consequence of change in M,, 



This image is compared with the reconstructed image of our observations, and the 
parameters of the PSF are varied until the best matching double-Gaussian is found. The 
final PSF of our LGS AO observations was: 0.17 and 0.80 arcsecs (FWHM) for the the 
narrow and broad Gaussian components and intensities of 0.53 and 0.47, respectively. 



Determining M, when the resolution is smaller than Hgph 

We observed NGC2549 with NIFS IFU centered on CO absorption features starting 
at 2.29 microns. We spatially binned the data using Voronoi method [9] and achieved 
the SN of ~ 60, while within the central arcsec the bins are no larger than 0.1 x 



0.1 arcsecs. We extracted stellar kinematics using the pPXF method [10]. Large scale 
kinematic data were obtained from previous observations with SAURON IFU [11] and 
were previously presented in [12]. We constructed orbit-based Schwarzschild dynamical 
models following the method presented in [13]. 

Fig. 1 shows grids of Schwarzschild dynamical models constrained using three dif- 
ferent kinematic data sets: combined SAURON and NIFS data, SAURON data only 
and NIFS data only. The best fit model using SAURON and NIFS kinematics gives 
M,=(1.4+j-3) X 10^ Msun- It is evident that when using only SAURON data it is not 
possible to determine the lower limit of the M,, simply because the resolution of the 
SAURON data is 1.7 arcsec, while Rsph-O.OS arcsec, using the above best fit model. Us- 
ing the NIFS dataset only it seems possible to determine even the lower limit, in spite of 

0. 17 arcsec resolution, but the uncertainty on M, is very large because, due to the small 
field-of-view, NIFS data alone are not able to constrain the total orbital distribution and 
give the right mass-to-light ratio. These are, however, well constrained using SAURON 
large field-of-view (FoV) observations. Combining the two data sets, the high resolution 
and large FoV, one can constrain Schwarzschild models and determine the M,. 

Fig. 2 illustrates that our NIFS observations, although at about 3 times lower resolu- 
tion than the estimated Rgph, are able to capture the change in the kinematics influenced 
by the SMBH. We illustrate this on the velocity dispersion: the a map of the best fit 
model is more similar to the observed data than a maps of models with a too small or a 
too big SMBH. To see this effect it was however, necessary to have high quality integral 
field data both at high resolution (NIFS) and with large FoV (SAURON). 
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